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GRAVIMETRIC INVESTIGATIONS ON THE NORTH AMERICAN DATUM 

(1972 - 1973) 

R. S. Mather 

Geodynamics Branch, Goddard Space Flight Center* 

SUMMARY 

All the available unclassified gravity data on the North American 
Datum (NAD) and in the surrounding oceans was assembled late in 
1972 for the investigation of the gravity field in North America and 
its relation to North American Datum 1927 (NAD 27), The gravity 
data in Canada and the United States was compiled on a common datum 
compatible with the International Gravity Standardization Network 1971 
(IGSN 71). The variation in the error of representation in the region 
is studied. Attempts are also made to study the correlation character- 
istics of gravity anomalies with elevation. 

A free air geoid (FAG 73) was computed from a combination of sur- 
face gravity data and Goddard Earth Model (GEM) 4 and this was used 
as the basis for the computation of the non-Stokesian contributions to the 
height anomaly. These non-Stokesian contributions as computed from 
the data sets available at present are found to occur with amplitudes 


*On leave of absence from tiie University of New South Wales, Sydney, Australia. 
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less than 3 m and with short wavelength in the Rocicy Mountains region. 
The resulting effects on determinations of the geocentric orientation 
parameters (geocentric datum shift) for NAD 27 are not found to be of 
significance at the 30 cm level. 

The geocentric orientation parameters obtained by this astro- 
gravimetric method are compared with those obtained by satellite tech- 
niques. The differences are found to be no greater than those between 
individual satellite solutions. The differences between the astro- 
gravimetric solution and satellite solutions GSFC 73 and GEM 6 are 
studied in detail with a view to obtaining a better understanding of these 
discrepancies. 
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GRAVIMETRIC INVESTIGATIONS ON THE NORTH AMERICAN DATUM 

(1972 - 1973) 

1. INTRODUCTION 
1.1 Preamble 

Studies made on the Australian Geodetic Datum (AGD) indicated that it ap- 
peared feasible to compile a world geodetic system (WGS) by comparing, in ef- 
fect, the elevations above ellipsoid as determined 

a. by gravimetric techniques above a geocentric ellipsoid; and 

b. by astro-geodetic methods above an ellipsoid defining the regional 
geodetic datum 

— y- 

for the evaluation of the geocentric orientation vector O (Mather 1970b). 

The experience obtained with investigations in Australia showed that com- 
parisons of the free air geoid determinations (Nf) obtained by the use of free 

air anomalies (Agj-) in Stokes’ integral 

R 

Nf = ;/f(^(/)Agfda (1) 

4tt'y 

(symbols and notation not otherwise explained in the text are described in Sec- 
tion 1.4) with astrogeodetic determinations on the regional geodetic datum, 
given by 

N. = -f f.lS 

.'geoid 


( 2 ) 
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where 


= Ik cosA^ 


(3) 


and 


Aj =* a; Aj - '/iTr-a (4) 

can be used to determine the geocentric orientation parameters A|j defining O 
in the local Laplacian triad at the origin of the regional geodetic system by 
(Mather 1970a, pp. 62 et seq, ) 

= A->h({j-{,) (5) 

where the elements hj^Alj^ of the vector h^AI^ at the regional geodetic origin 

( 0 Q> related to the equivalent parameters hj (|gj - |^j) defining the vec- 

•> 

tor h(|g - Ig) at the j-th station (^,X) on the regional datum by the elements in 
the matrix 


A = 


sin 0Q 


AX beii^ given by 


sin 0 sin AX 


- sin sin AX cos AX 


cos 0Q sin 4> - 
sin cos tp cos AX 


- cos 0 sin AX 


AX = Xq - X 


and 


sin cos 0 - 
cos 0Q sin 0 cos AX 

- cos 0Q sin AX 

sin sin 0 + 
cos 0Q cos 0 cos AX 


( 6 ) 


hj = -(p + h); h, = -(v + h): II 3 = 1 (7) 

all other notation being described in Section 1,4, Strictly compatible compari- 
sons of geocentric quantities defined from gravimetric considerations (denoted 
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by the subscript g) with equivalent quantities on the regional datum obtained 
from astro-geodesy (subscript j, ) are the following; 

?at = ^al -O.17(sec)h(km)sin20+A|j 

^g2 ~ ^f2 ^c2 ’ ^a2 “ ^a2 ^^2 

( 8 ) 

Ig3 (~ b) - |jj 3 (“ h) h |3 + ^ dz- ^ 

geoid geoid 

where are the two components of the astro-geodetic deflection of the vertical 
on the regional datum, hp is the height of the levelling datum { bove the regional 
reference ellipsoid, the component of the astro-geodetic deflection of the 
vertical in the line of levelling of length dfi along which the difference of ortho- 
metric elevation is dz. It should also be noted that exact equivalence is obtained 
in the third equation at 8 only when dz and refer to the same levelling route. 

On ignoring these relative subtleties in the definition of A^ 3 , a solution was 
made for the geocentric orientation parameters A^j for the AGD from a direct 
comparison of the free air geoid with an astro-geodetie determination based on 
Equation 2 (Mather 1970a) as it was assessed that the indirect effect for a region 
of limited topographic variation like Australia was less than 60 cm (Fryer 1970). 

It was also argued that tliis procedure could be the basis for the assembly 
of a world geodetic system (WGS) without (eventually) having recourse to a satel- 
lite model of the Earth's gravity field (Mather 1972a). Definite conclusions on 
the effective practical utilization of such a system on the basis of the Australian 


study alone were limited by 
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i, the extent of the AGD (restricted to only 1-1/2% of the global siu’face 
area); and 

ii, the relatively flat terrain ove.’ wliich the magnitude of the non- 
Stolcesian term tended to be zero. 

The study was therefore extended to the North American Datum (NAD) 
which covered about 2-1/2 times the surface area of the AGD, with a view to 
investigating 

a. the significance of the non-Stokesian contributions to determinations 
of the geocentric orientation vector; and 

b, the stability of appropriate harmonics as determined from "surface 
fitting," if the astro-geodetic data on NAD was sufficiently dense to 
warrant such studies. 

1.2 Data Available on the North American Datum in 1972 

The basic sources of information for the assembly of a gravity data banlt 
on NAD were the following: 

• The holdings of the United States N’ttional Geodetic Survey, based 
largely on unclassified material assembled by the Defense Mapping 


Agency Aerospace Center. 
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• Data made available by the Gravity Division, Canadian Department of 
Energy, Mines and Resources, Ottawa. 

1‘his information was supplemented by gravity data collected in the sur- 
rounding oceans on various surveys and obtained for the Geodynamics Branch, 
Goddard Space Flight Center, by Computer Sciences Corporation, Silver Spring, 
Maryland. Other gravity data for these oceanic areas were in the form of 1° x 
1“ area mean free air anomalies (ACIC 1971; Taiwan! et al. , 1972). The basic 
gravity data was made compatible with International Graviiy Standardization 
Network 1971 (IGSN 71). This required that a correction of +2.0 mGal be made 
to the Canadian data to bring it to a common datum with the U.S. gravity data 
{Tanner , 1972). 

The latest low degree surface spherical harmonic model available in the 
second half of 1972 for the Earth's disturbing potential was that prepared at 
Goddard Space Flight Center— Goddard Earth Model (GEM) 4— from the combina- 
tion of data obtained from orbital analysis with surface gravity information. 

This model was used along with surface gravity data from the sources mentioned 
above to provide a continuous representation of the global gravity anomaly field 
as described in Section 3.1. 

The effect of the terrain undulations on solutions of the geodetic boundary 
value problem arise entirely in the non-Stokesian term. An elevation data banit 
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is required for the meaningful evaluation of these effects. The elevation data 
available for this study in 1972 were parts of the following data sets: 

• 5’ X 5' mean elevations for the NAD region below parallel 65 °N and 
made available by the Defense Mapping Agency Aerospace Center 
(Czai’necld, 1970); and 

• 1° X 1° mean elevation estimates for the world prepared at the 
University of California (Kaula et al. , 1966). 

1.3 Computational Considerations 

The computations proceeded on the following lines: 

i. The available surface gravity and elevation data was processed to 
obtain e.itimates of the area mean free air anomalies and elevations 
on a tenth degree equi-angiUar grid covering North America. 

ii. Prediction techniques were used to estimate values of the free air 
anomalies in unsurveyed regions covered by this tenth degree grid 
in locations where the surface gravity information warranted this 
approach. For further details, see Section 3.1. 

iii. The resulting continuous data set was used to compute consistent 
1/2°, 1° and 5° equi-angular area means for all regions of relevance 
as indicated in Table 1, the surface data being augmented by GEM 4 


model values in distant areas and wherever warranted. 
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iv. This data set was used to compute the free air geoid (Nf) for North 
America (FAG 73) using the general techniques previously applied to 
the AGD as described in (Mather, 1970a). 


V. FAG 73 was used along with the tenth degree equi-angular elevation 
mean values to compute the non-Stokesian contribution to the height 
anomaly f given by (Ibid, p. 86) 


r = Nf+N^ 

where 

Nf = 1 ( (W„ - U„) - R M{Agf } ) + ^ ;; «*) Agf da 

and 


R2 


1 /R 


Ng = ff — — sin i// 


dh hp-h\T 


dr 


R 


tan/3„ 


da(+o{fblJ if'/ 2 (h -h)2/r^<f) 


where 


— = tan cos A'' 

A'/ = A„; A'' = 'Att-A, 

all other quantities being described in Section 1.4. 


( 9 ) 


( 10 ) 


( 11 ) 


(12) 

(13) 


The first term on the right in Equation 10 is of relevance only if 
Agj- is a global sample measured at the surface of the Earth. This 


term is ignored in the present study as this is not currently the case. 
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The condition under which the relation for N^. , given as Equation 11, 
is applicable can be said to be satisfied when the elevation data input 
is in the foi’m of a tenth degree area mean data banic. As discussed 
in (Mather, 1972a), the second term on the right of Equation 11 can 
contribute up to 10% of the Stokesian effect in those momitainous con- 
tinental regions where the sign of the deflections of the vertical is 
strongly correlated with that of the topography. The magnitude drops 
off to "lero in predominantly oceanic areas. Its total magnitude mider 
xmfavourable conditions is not ejqDected to contribute in excess of the 
order of tlie flattening to the magnitude of f for i//>5° even though a 
tendency for negative values of tan appears likely over the 2-3% 
of the Earth's surface area where the great momitain ranges occur. 
This is due to the tendency for the wavelengths of such contributions 
to be rather small. 

The first term on the right of Equation 11 is a function of r“^ and 
therefore damps out rapidly (i.e. , < o{f}) for l^'> 1-1/2° . Its magni- 
tude is dominant in mountainous coimtry and the sign of the contribu- 
tion is variable, depending on the relation of the surrounding terrain 
to the elevation of the point of computation. The apparent instability 
of this term as is offset by the fact that sin i//(dh/dr) -*■ (hp - h)/R 
as This enables Equation 11 to be evaluated from the limited 

10 km mean elevation data banic available without seriously undermining 


the stability of computations. 
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It is obvious that T has to be approximated by the free air geoid Nf 
and the components of the deflections of the vertical ) are computed 
by using free air anomalies in the Vening Meinesz integrals 

Further Information on the techniques of computation is given in 
Section 3. 

vi. Geocentric orientation parameters A^. are computed for NAD 27 both 
in the Laplacian triad at the Meades Ranch Origin as well as in rela- 
tion to a three dimensional Cartesian coordinate system Xj described 
in Section 1.4. Numerical values for A^j are obtained by comparing 
either Nj. or f as obtained above with asti'o-geodetic values of N^ as 
defined in Equation 2, using Equations 5 through 3. The results ob- 
tained are discussed in Section 4. 

1.4 A Guide to Notation 

A = parameter associated with azimuth (Equations 4 and 13) 
a = equatorial radius of reference ellipsoid 

dfi = distance between terminal bench marks in loops of levelling 
(Equations 2 and 8) 

dz = increment in orthoraetric elevation between terminal benchmarks 
in loops of levelling 

d (7 = element of surface area on unit sphere 
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•f 

m 


G 


g 

h 



m 


Nc 

Nf 

-»• 

O 



r 


r 


o 


= error of representation for a n° x n° equi-angular element of sur- 
face area 

= flattening of reference ellipsoid 

= Stokes' function = 1 + cosec Vi\j/ - 5 cos ij^ - 6 sin '/up - 3 cos 

(15) 


p log (sin Yip + siir Vip) 

- gravitational constant 
= observed gravity at the surface of the Earth 

= elevation above ellipsoid 

AW / , AW 

= normal height I 1 - (1 + 1 + m - 21 sin- p) — 

T V 



+ o{f3} 

(16) 


= global mean value of X 

= aco^ /7g 

= astro-geodetic "geoid" 

= non-Stokesian contribution to the height anomaly (Equation 11) 
= free air geoid (Equation 10) 

= geocentric orientation vector defining the displacement of the 


origin of the regional geodetic datum from the geocentre (Earth's 
centre of mass) 

= order of magnitude of the largest terms not considered is that of X 
= mean radius of the Earth 

= distance of the element of surface area dS from the point of com- 
putation P 
= 2R sin Vip 


(17) 



disturbing potential 

potential of reference system on the equipotential reference 
ellipsoid 

potential of the geoid (unlmown) 

geocentric three dimensional Cartesian coordinate system with 
axis passing through CIO pole and the Xj X3 plane defining tliat of 
zero longitude 

azimuth from the point of computation 

groimd slope; use with subscripts denote north and east components 
global mean value of normal gravity 

equatorial normal gi’avity ' 

value of normal gravity on equipotential reference ellipsoid 

AW / AW \ 

gravity anomaly = g-j - 2 — (l+f+m+ -2fsin^ 0 + o{f2] (18) 

“ a \ 2 a 7 j 

free air anomaly = g -7^ + h^"’^ (19) 


c 

geopotential difference with respect to the geoid = - I S dz (20) 

•'geoid 

components of the geocentric orientation vector O in the Laplacian 
triad at the regional geodetic origin (Equation 26) 
height anomaly 

component of the deflection of the vertical in azimuth a 
longitude positive east 

radius of curvature of reference ellipsoid in prime vertical normal 


section 
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= conipoiiont of doflQction of the vertical In the meridian (k =1) and 
prime vertical (k= 2) directions, positive if outward vertical is 
nortli or east of normal 

p = radius of curvature of reference ellipsoid in mex’idian normal 
section 

(j) - latitude, positive north 

^ - angular distance between tlie variable element of surface area do 

and the point of computation 
CO = angular velocity of rotation of tlie Earth 

Significance of Subscripts 

„ = astro-geodetic values referred to the regional geodetic datum 

= non-Stokesinn contribution for conversion of free air geoid value 
to physical surface/toliuroid system 
p = free air geoid 

g = geocentric (gravimetidcally determined) values 

= values at the regional geodetic origin 
p - values at the point of computation 

o = values at tiie element of surface area do 

Note ; 

1. Repeated subscripts in a product indicate summation over all possible 


values. 
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ii. Repeated subscripts on both sides of an inequality indicate as many 
equations as possible values of subscript. 

iii. Use of roman letters for variable subsci’ipt indicate three values; the 
number of possible values when greek lettei'‘s are used is two . 

2. GRAVITY DATA PROCESSING 

2. 1 The Area 

All the available gravity data in the area 

0° <0<9O°; -180° 

was processed along the lines described in tlois section. In the first instance, 
the area was divided into 81 10° x 20° regions. The continental area of the 
North American Datum (NAD) was covered by the twenty nine regions shown in 
Figure 1. The gravity data in each of these regions was processed separately. 
Additional input into the basic processiiag routines was the equivalent equi- 
angular mean elevation data banic for North America described in Sections 1.2 
and 1.3. 

After all the available gravity data had been sorted and converted into 
free air anomalies consistent with IGSN 71 in the appropriate tenth degree equi- 
angular squares, area mean values of the free air anomaly (Ag^) and gravity 
station elevation (h^ ) were computed for each square from the available sample. 
The data was then adjusted to the mean elevation (h) of the tenth degree square 
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obtained from the continuous elevation data set referred to above, according to 
the relation 

Agf = Ag,-0.112("»°='l)(h(m)-h(nO) (21) 

It was believed that the adoption of such a procedure should eliminate any tendency 
for gravity data samples to be biased as a consequence of data being collected in 
the more accessible valleys. For further comment on the validity of this as- 
sumption, see Section 2.2. 

The distribution of gravity data on the North American continent which was 
available for the present set of calculations is shown in Figure 2. The represen- 
tation of the free air anomaly data set as contoured from one degree equi -angular 
area means is shown in Figure 3. The region is one of predominantly negative 
free air anomalies with positive values occurring primarily in the elevated re- 
gions. As mentioned earlier, all the Canadian data is 2mGal larger than in pre- 
vious representations which were not on IGSN 71, all data having been adjusted 
in accordance with Equation 21 in tenth degree equi-angnlar squares. 

2.2 Characteristics of the Free Air Anomaly Data 

Noi’th America is a region of varying topography with the Rocky Momitains 
dominating the western part as illustrated in Figure 4. This provided the oppor- 
tunity for a closer look at the characteristics of free ail' anomalies in relation 
to topography. A histogram showing the distribution of one degree area mean 
elevations for equi-angular squares as a function of elevation is shown in Figure 5. 




It is the writer's opinion that the pre-requisite for a successful evaJuation 
of the height anomaly by quadratures methods from irregularly distributed data 
is an appreciation of the variability of the gravity anomaly field. A study of 
Figure 3 indicates that the free air anomaly field is questionably holonomic in 
two dimensions at the surface of the Earth. It also points to some correlation 
of the magnitude of free air anomalies with elevation though this tendency could 
well be masked by other factors. In the present investigation, the gravit7 data 
was sorted initially into tenth degree equi-angular squares. Area mean free air 
anomalies for larger half degree, one degree and five degree equi-angular 
squares were computed from the basic tenth degree values after prediction of 
estimated free air anomalies in unsurveyed areas using techniques described in 
(Mather, 1970a, pp. 65 et seq.). In tliis manner, all data sets were fully repre- 
sented and were inter-compatible. 

The variability of free air anomalies in these various squares can be 
quantified by the error of representation E{Ag}j^^o for a m° x m“ square. 

/ M N; M Nj 

2 (Agij-Agj)2 cos0y/ 2] 2 COS0JJ (22) 

= ) i=ij = l / 

where Nj is the number of readings in the i-th m° x m° square and M is the 

total number of such squares. 

It was initially hoped to evaluate definitive values of E{Ag}j^o for the 
square sizes shown in Table 1 for the North American region. In practice, how- 


ever, it was difficult to achieve tliis goal for square sizes smaller than 0. 5° due 
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to the irregular sample distribution and the general paucity of data at intervals 
smaller than 10 km, hi the case of squares where data was available with a 
greater density, irregular concentrations resulted in abnormally low values 
for E{Ag}jj jO and E{Ag}p qs® ■ These results are therefore not presented. 

E{Ag}^o as obtained for North America from the data sample available 
for the present study is given in Table 2 for m® =0.5, 1 and 5. hi view of the 
variation in the surface area of equi-angular squares with latitude, a plot of 
E{Ag}g jO against latitude is illustrated in Figure 6, together with a plot of the 
mean gravity station elevations and the mean elevations for all the half degree 
squares along the stated parallel of latitude. There appears to be no significant 
variation in the magnitude of E{Ag}g with latitude except in the south where 
the proportion of mountainous terrain is greater as indicated by the higher values 
for the mean elevation of the topography. This is principally due to the fact that 
the maximum possible distance between readings in a single surface element 
always remains the same (0.5° in latitude). 

A more uignificant effect is the possibility of correlation of E{Ag}|^»with 
elevation as illistrated by the values for E{Ag}j» in Table 2. A histogram show- 
ing the frequency of occurrence of the various numerical values of E{Ag}, “ 
in North America is illustrated in Figure 7. TMs figure illustrates that a wide 
range of values exists for E{Ag}jO on NAD with a modal value of around ±10mGal. 
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The information illustrated in Figure 6 is made more meaningful if the two 
factors mentioned above are taken into account. The much less significant lati- 
tude effect is easily accounted for by using Equation 22. The more dominant 
effect of elevation on E{Ag}j^o is studied by classifying the correlations that ex- 
ist between such values and elevation. Results for m = 1 are illustrated in Fig- 
ure 8. The data analyzed was restricted to those one degree equi-angular 
squares containing more than 40 tenth degree equi-angular area means based on 
observations (i.e. , 40% of the total possible representation). 

It can be seen that no meaningful trend indicating a correlation of the mag- 
nitude of E{Ag}jofor free air anomalies with elevation can be obtained due to 
three reasons: 

i. There are significantly more squares with lower elevations than higher 
ones (Fig. 5) 

ii. E{Ag} is correlated not so much with elevation as with ruggedness of 
topography. (Thus E{Ag}for a high plateau should be no more suscepti- 
ble to larger than average values than is the case for a low plain. ) 

lii. Values of E{Ag} for coastal areas tend to be disturbingly large, as il- 
lustrated in Figure 9. Tliis could be due to the greater variability of 
terrain/sea bed topography in continental margin areas close to the point 
of gravity measurement. It is hoped that this feature is not due to tlie 
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fact that oceanic data is usually supplied by sources different from 
those providing land data and/or associated incompatibilities in the 
different data acquisition systems. 

Despite the limited samples for the less commonly occurring elevated areas 
In the NAD region, it is nevertheless useful to study the values of E{Ag}(j 5 ° 
as classified by elevation and shown in Table 3. The value for the elevation 
range 0 ^h<1000 is similar to that obtained for Australia (Mather, 1967, p. 131), 
In assessing the overall pattern of variation of E{Ag}p for North America as 
categorized in Table 3, it would not be unreasonable to expect the variability of 
the free air anomaly field to be relatively smaller in coastal plains and deep 
oceans, with E{Ag} talcing greater magnitudes in the vicinity of continental shelf 
and in mountainous regions. Figures given in Tables 2 and 3 suggest that values 
of E{Ag} in the vicinity of shelf areas and in elevated regions could be larger 
than the average value for coastal plains by a factor of between two and three. 

Such considerations have an important bearing on the procedure to be 
adopted in the prediction of values to represent unsurveyed areas when seelcing 
a model for the fine structure of the free air anomaly field. It is common 
practice to eliminate height correlation effects by converting free air anomalies 
to Bouguer anomalies prior to prediction (e.g, , Molodenslcii et al. , 1962, p. 

179; Heiskanen ^ Moritz, 1967, p. 281; Mather, 1967, p. 131). An attempt 
was therefore made to study the height correlationcharacteristics of the available 



free air anomaly data on NAD by linear regx’ession analysis of the data in the 
three larger square sizes {i.e. , half degree, one degree and five degree). 

In the case of five degree equi-angular squares, the results indicated that 
the gradient of a plot of free air anomalies against elevation varied consider- 
ably and seldom approached the expected value of O.llmGal m“i . While the 
tendency towards correlation of free air anomalies with elevation increased in 
the case of data confined to a one degree equi-angular square, again, the gradi- 
ent did not approach the value quoted above with any degree of certainty. On the 
basis of the data studied, it must be concluded that if there were a positive cor- 
relation of free air anomalies with elevation, it is likely to be masked by other 
factors when the distance involved exceeds 50 km (the minimum distance con- 
sidered in the present investigation). It can therefore be concluded that any 
tendency for free air anomalies to be correlated with elevation on a regular 
basis capable of analytical representation, has not been established from the 
available data sample on the NAD in the case of distances in excess of 50 Ian. 

Due to the time constraints on this investigation, the statistical analysis 
of the free air anomhly data was done after the assembly of the data set de- 
scribed in Section 3,1. Consequently, the prediction of values to represent 
unsurveyed areas was based on the conventional assuvaptions referred to in the 
previous paragraph plus one, as controlled using the tecimique described in 
(Mather, 1970a), As a result, large positive anomalies tend to occur when 
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predicting in mountainous regions. It is unlikely that tliis affects the results 
outside the local area in excess of ±30 cm in view of the short wavelength of such 
biases (Mather, 1974, p. 103). 

The writer plans to implement prediction procedures in future invescigations 
on the basis of findings in tliis study. Irrespective of the exact prediction pro- 
cedure adopted, they can be broadly classified as follows: 

i. Choose the basic surface area size in which the gravity anomaly data 
is to be processed (as a percentage of the total surface area of the 
Earth). 

ii. Establish a variability characteristic (e.g. , error of representation) 
for this block size in the case of both gravity anomalies and elevations 

a. for each individual area; and 

b. for the entire sample of such areas. 

Tliis characterization of the variational behaviour of the free air anomaly field 
in each basic area can be embodied in either a statistical measure such as co- 
variance (about zero mean for each area at i.) or the number Uj^a^ 
monic coefficients to be used in the analysis of the available data. In the latter 
case, the value of would also be influenced by the number of readings 
available in each basic area, ha general, the larger the error of representation, 
the greater the value of for a given data distribution. For further 
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discussion on the relation between and the distribution of data in a 

region, see (Mather, 1967, p, 133). 

3. COMPUTATIONS 

3.1 The Free Air Geoid (FAG 73) 

The gravity data referred to in Section 1.2 was used to define a continuous 
free air anomaly field on the following basis: 

In regions where point data was available within a 1° x 1° square 

Prediction was preceded by the conversion of all free air anomaly means for 
tenth degree equi-angular areas to Bouguer anomaly equivalents using the tenth 
degree elevation data banlc. Two dimensional harmonic series were fitted to the 
residuals from the GEM 4 model free air anomalies using the procedure de- 
scribed in (Mather, 1970a, pp. 65 et seq.). This type of approach was used in 
the twenty nine regions shown in Figure 1 despite the reservations expressed in 
Section 2.2 about assuming predominantly height correlated variations of free 
air anomalies over distances in excess of 50 km. The final set of free air an- 
omalies was generated after the completion of prediction by reversing the Bouguer 
correction and restoring the GEM 4 model. 

In regions where point data was not available within a 1° x 1° square 

These regions were largely confined to the surrounding oceanic areas where 
the required field representation in terms of Table 1 were either in the form of 
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ha).f, one or five degree equi-angular area mean free air anomalies. The follow- 
ing hierarchical structure was adopted for this type of region: 

1. Lamont 1° x 1° area means (Taiwan! et al. , 1972), 

2. DMAAC 1“ X 1° area means (ACIC, 1971). 

3. GEM 4 model anomalies (Lerch et al. , 1972), 

A consistent and continuous representation of the free air anomaly gravity 
field was constructed on this basis for the ai'ea defined in Section 2.1. Thus, 
the tenth degree equi-angular data set extended only up to 3“ from the coastline, 
while the half degree data set extended up to about 8° from the coast. All data 
was maintained on Geodetic Reference System 1967 - GRS 67 (LAG, 1971). The 
fact that tiiis system differs from current estimates of best fitting Earth ellip- 
soid parameters (e.g. , Lerch et al. , 1974) was not considered to be of signifi- 
cance for the following reasons: 

a. Zero degree terms were not considered in the gravimetric solution for 
reasons given in Section 1.3(v). Thus the exact value of a and GM used 
have only marginal effects on the solution smaller than the noise pro- 
duced by the errors in the input data. 

b. The only reference model parameter which has an appreciable non-7.ero 
degree effect is the flattening f. As the best fitting value at the present 
time is held to be around 1/298.255 (e.g. , Ibid) instead of the 1/298.247 
implicit in GRS 67, the resulting effect is of order 10“'* which can be 
ignored in solutions seeking a precision less than ±5 cm. 
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While the free air anomalies based on point values were effectively refer- 
enc^d to IGSN 71, oth^ data hi the.fprnj onebegre-e,^^^ were implic- 

itly corrected to this system by applying the recommended correction to the 
Potsdam datum (lAG 1971). The resulting free air geoid for Nortli America is 
called FAG 73 in the text and was obtained by tlie use of Equation 1 and the area 
sub-divisions specified in Table 1. A contour representation of FAG 73 based 
on a sample plot of values on a two degree grid is illustrated in Figure 10. 

3.2 Comments on the Precision of FAG 73 

The computational procedures used in the preparation of FAG 73 were de- 
signed in order that the precision attainable from an error free data set was 
±30 cm. The result illustrated in Figure 10 cannot be expected to have achieved 
this precision for the following reasons: 

• Any errors in GEM 4 which was effectively used to represent the dis- 
tant zones, would be reflected in the results. This type of error is 
tentatively estimated as slowly varying at the 1-2 m level. 

• The oceanic repx’esentation, especially in the Pacific, is not based on 
any systematic examination of the data by the investigator, the num- 
bers given in the data sets descx’ibed in the second paragraph of Sec- 
tion 1.2 being accepted at face value. Nevertheless, both these data 
sets are in reasonable agreement with the data set prepared by the 
writer in areas of overlap on allowing for the different methods used 
in the preparation of the latter. 
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• Prediction errors in scantily surveyed areas (see Pig. 2) could also 
degrade the quality of the result. 

Errors in excess of the norm could tlierefore be expected for determinations in 
Alaska, the west coast of North America and Mexico. This expectation appears 
to be substantiated when comparing FAG 73 with astro-geodetic determinations 
after allowance for datum translation. Such a comparison is illustrated in Fig- 
lu'e 11. Details of how such comparisons are effected could be ascei’tained by 
reference to Section 4. 

In assessing the nature of the discrepancies between gravimetric and astro- 
geodetic solutions as illustrated in Figure 11, it should be borne in mind that 
FAG 73 and the asti’O-geodetic imdulations detenuined by the application of 
Equation 2 are not sitrictly comparable for the reasons set out in Equation 8. 

The extent of the discrepancy is likely to be exacerbated in elevated regions and 
this is borne out by the largest non-coastal discrepancy between the solutions 
occurring in the most elevated region on the datum in Colorado. About 1-1/2 m 
of the discrepancy is accounted for by the non-Stokesian contribution to the height 
anomaly (see Fig. 12). Another factor contributing to this discrepancy is the ap- 
proximation of orthometric elevation in the area by the sum of the orthometric 
height differences as observed along the line of levelling. The resulting effect 
in an area with elevations in excess of 2000 m could be as large as 2 m. 
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The largest discrepancies occur along the north west coast of North 
America. The quality of neither the graviiy data (Fig. 2) nor the astro-geodetic 
data (Fig. 14) in this area warrants further speculation on the underlying causes 
of this discrepancy. 

FAG 73 shows a less accentuated hump over the Rocl^y Mountains in Colorado 
that the Vincent-Marsh determination (Vincent & Marsh, 1973) or— as is more 
relevant— the GEM 4 geoid (Lerch et al. , 1972). This is attributed to the effect 
of the relatively denser representation of the predominantly negative free air 
anomaly field afforded by the data set used. The other prominent feature is the 
geoidal low over Hudson Bay wliich is usually shown as being lower than 50 m. 

The low illustrated is at least two meters higher due to the +2 mGal adjustment 
made to all the Canadian data as described in Section 1.2. 

For statistics on tests between astro-geodetic solutions and FAG 73, see 
Table 4 and Section 4. 

3 . 3 Non-Stokesian Contribution to the Height Anomaly 

The non-Stokesian contribution to the height anomaly as estimated for North 
America is illustrated in Figure 12. As explained in Section 1.3(v), the magni- 
tude of Np was computed using as input 

• FAG 73 computed using Equations 1 and 14; and 

• the tenth degree equi-angular mean elevation data bank 
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in Equation 11, As expected, the largest contributions occur in the Roclty 
Mountains area, never exceeding ±3m in magnitude. The unexpected feature 


was the occurrence of both positive and negative magnitudes, dominant magni- 
tudes of the latter occurring on the oceanic flanlcs of the Roclcy Moimtains. As 
these non-Stokesian contributions have iiinited amplitude and occur with relatively 
short wavelength, their effect on determinations of the geocentric orientation 
parameters are fomid to be negligible. TMs is discussed at length in Section 4, 

4. COMPUTATION OF THE GEOCENTRIC ORIENTATION VECTOR FOR 
THE NORTH AMERICAN DATUM 1927 

4.1 Transformation of Astro-Geodetic Data From Clarke 1866 to Reference 
Ellipsoid 1967 

The use of Equations 5 to 8 for the definition of the geocentric orientation 

parameters A^j in the local Laplacian triad at the Meades Ranch origin of NAD 

— ^ 

27 and hence the geocentric orientation vector O is possible only if the astro- 
geodetic data were referenced to Reference Ellipsoid 1967 (a^ - 6 378 160 m; 
f"’ = 298.247) placed concentric with the Clarke 1866 ellipsoid (a = 6 378 206.4m; 
f“’ = 294. 979) in Earth space. This is achieved by allowing for the corrections 
(e.g. , Mather, 1968, p. 292) 

50 = sin 20 ^df+ f ^ -df ^ + fdf cos- 0 + o{i'»}j (23) 

to 0 and 


5N = da + a(df sin^ 0-50 tan 0) + o{f^a} 


(24) 
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to the height h above the ellipsoid, where 

da = a^ - a and df = - f (25) 

The geometry of the system is illustrated in Figure 13 which shows the Clarke 
1866 and the 1967 Reference Ellipsoids defining NAD 27 in relation to GRS 67. 

4.2 Geocentric Orientation Parameters for NAD 27 from Astro-gravimetric 
Comparisons 

4,2.1 Astro-Geodetic Data 

Astro-geodetic data on the NAD is available in one of two forms: 

i. Astro-geodetic deflections of the vertical distributed as shown in 
FigTire 14, 

ii. Astro-geodetic "geoids” hereafter referred to as astro -geodetic un- 
dulation^^ for want of a more accurate description and 

computed using equations of the type at 2 and 4. Solutions of this t 3 rpe 
were produced by 

a. the U.S. National Geodetic Survey for selected loops of astro- 
geodetic stations in the United States (Rice, 1972); 

b. the Canadian Geodetic Survey for south-eastern Canada (Ney, 1952); 
and 

c. the U.S. Array Map Service as a contour map for the whole of North 
America (Fischer et al, , 1967). 
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4.2.2 Comparison of Gravimetric and Astro-Geodetic Determinations 

Geocentric orientation parameters obtained from the comparison of gravi- 
metric and astro-geodetic determinations of the separation vector can be de- 
scribed as being largely achieved by the process of surface fitting and have al- 
ready been used to provide a geocentric orientation for the Australian Geodetic 
Datum (AGD). These comparisons also provide a means of assessing the pre- 
cision of both the gravimetric and astro-geodetic determinations (Mather, 1970a; 
Mather, 1972b). 

Geocentric orientation parameters for the major geodetic datums are freely 
available in the form of components on a geocentric Cartesian coordinate sys- 
tem. They are obtained by comparing the coordinates of satellite tracldng sta- 
tions in geocentric solutions with equivalent values assigned to these same sta- 
tions as part of regional geodetic surveys on the local datum. Results for NAD 27 
are given by several investigators (e.g. , Anderle, 1974; Gaposchldn, 1973; 

Lerch et al. , 1974; Marsh et al., 1973; Merry & Vanicek, 1974; Mueller, 1974; 
Schmid, 1974). The technique used in this present study is different from satel- 
lite based determinations in several important respects: 

a. The geocentric orientation parameters are defined by fitting estimates 
of the shape of the same surface as determined by different methods— 
one in relation to the geocentre and the other in relation to NAD 27— on 
allowing for the difference in shapes of the two reference surfaces (see 


Section 4.1). 
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b. The determination is based on comparisons made over the entire extent 
of the datum rather than at a few selected tracldng stations whose geo- 
metrical distribution across the datum could well be irregular. 

c. The results obtained from satellite solutions and summarized in Tables, 
are obtained by comparing coordinates of the same points on the two 
different datums on the basis of a seven parameter fit (3 translational - 
AXj, 3 rotational - Wj and one for scale -n ). The astro-gravimetric 
method, on the other hand, defines directly a three parameter trans- 
lational shift as explained in the Appendix. 

d. Zero degree effects, if any, in the height anomaly are not included in 
the solution. 

e. The size of the ellipsoid and GM are not critical factors in the solution, 
provided the differential geometry is correctly allowed for. 

Several solutions were made for the geocentric orientation parameters Afj 
for NAD 27 by comparing gravimetricaily determined height anomalies or co- 
geoid heights for FAG 73 on a geocentric ellipsoid witli astro-geodetic undula- 
tions on NAD 27. Unfortunately, for the very reason that a re-definition is 
currently being undertaken for the NAD (see Canadian Surveyor 28(5), 1974, 
for details) the astro-geodetic data available in 1973 could be expected to be 
subject to limitations when used in surface shape comparison methods for the 
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definition of tlie geocentric orientation vector. It has been noted that errors as 
large as 15ppm could exist in the Canadian network (Lilly, 1960). The existence 
of such uncertainties would almost certainly limit the precision which could be 
achieved by the use of astro-gravimetric methods at the present time. 

Nevertheless, it was considered worthwhile to evaluate the geocentric 
orientation vector for NAD 27 for the following reasons: 

• The precision of the gravimetric determinations could be estimated by 
comparison with astro-geodetic determinations after translation of 
datums. 

• The effect of the non-Stokesian contribution to the height anomaly on the 
determination of the geocentric orientation vector could be evaluated. 

• It is of interest to define the geocentric orientation vector using a 
technique not sensitive to scale apart from zero degree contributions 
to the height anomaly. It is not expected that the latter would exceed 
3-6 m (Mather, 1970b, p. 98). Any such effect would, however, only 
influence the determination in the radial direction at Meades Ranch . 

• Differences between the scale of satellite solutions and that of the re- 
gional geodetic network not dependent on the choice of ellipsoid will, 
however, be reflected in the results. 
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4.2.3 Solutions Using the Free Air Geoid (FAG 73) 

The different types of comparisons between astro-geodetic and gravimetric 
solutions made for geocentric orientation parameters are listed in Table 4. 
Comparisons were effected between gravimetric determinations and two types 
of asti’o-geodetic undulation repi-esentations: 

1. the astro-geodetic undulation map produced by the U.S. Army Map Ser- 
vice (Fisher et al. , 1967), The version used was digitized by inter- 
polation on a one degree equi-angular grid and called AMS 67. 

2. In this representation, the astro-geodetic data was sampled on a one 
degree equi-angular grid. These grid corners were represented 
wherever possible by the nearest astro-geodetic station. Approxi- 
mately 800 of the available astro-geodetic stations on the NAD were 
used in this set on this basis. Undulation values at 364 of these sta- 
tions were obtained either from the astro-geodetic levelling results 
oftheU.S. National Geodetic Sui'vey (NGS 72) prepared by Rice and 
co-workers (Rice, 1972) or the results for south-eastern Canada re- 
ported by Ney (1952). (Note; There are over 3000 deflection stations 
available in North America but only about 800 were selected to provide 
an "equal area" coverage of the continent.) 

Solutions based solely on undulation comparisons are classified as Class A 


in Table 4, The geocentric orientation parameters A^» Atj and AN in the Laplacian 
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triad at Meades Ranch are related to Equation 5 by 

A^j = A| = -50: A^2 = At? = -SXcoscj)^-, A^j = AN (26) 

where 60, 5X and AN are the changes in latitude, longitude and normal displace- 
ment at the geodetic origin of the NAD required to give equivalent geodetic co- 
ordinates on a geocentric ellipsoid. It should, however, be recognized that 
comparisons between the gravimetrically determined undulations on the one hand 
and equivalent astro-geodetic quantities on the other, should be preceded by a 
transformation of the latter from values referred to the Clarke 1966 ellipsoid 
to corresponding values on the concentric ellipsoid with dimensions equivalent 
to those used in Geodetic Reference System (GRS) 1967. Details are given in 
Section 4. 1. The values 

A^q ” -7.5 arc sec and AN^ = 48m (27) 

should be subtracted from values in columns 5 and 11 in Table 4 to obtain equiva- 
lent values for shifting NAD 27. 

Class A solutions have been effected on the following basis in the case of 
comparisons between undulations only of 

• PAG 73 on the one hand and AMS 67 (representation on a 1° x 1° grid) 
on the other (Table 4, Rows 1 & 7); 

• PAG 73 and astro -geodetic undulations from solution NGS 72 for USA 
as supplemented where relevant by Ney's solutions for south-eastern 
Canada at selected astro-geodetic stations as described above (Table 


4, Rows 3 & 6). 
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The nature of the data set AMS 67 is different fi’om that of NGS 72 and Ney. The 
values comprising the former are obtained by interpolation from a contour map 
wliile the latter solutions are based on values at points in loops of astro-geodetic 
levelling which are therefo'’e more reliable. This is borne out by comparing the 
values in column 12, Rows 6 & 7 of Table 4, wMch are based on comparisons 
over substantially the same area. 

The geocentric orientation parameters obtained by tliis technique of astro- 
gravimetric comparisons are, as expected, a function of the area over which 
the comparisons are made. The parameter least affected is AN which has an 
overall variation of less than Im for all Class A solutions, irrespective of the 
extent of the area over which comparisons are made. 

On application of the geocentric orientation parameters obtained from the 
use of Equations 5 through 8, to transform the astro-geodetic determination to 
GRS 67, it is possible to compare values of 

• the surface undulation; and 

• the deflections of the vertical 

as obtained by these two independent methods to provide a measure of the suc- 
cess with which the orientation has been achieved. 

The goodness of fit can be characterized by the root mean square (rms) 
residuals in the comparison of the undulations (Table 4, Column 13), 

(Table 4, Colunm 6) and (Table 4, Column 9) in the deflections of the vertical. 
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It wovild appear at first glance that over similar areas (e.g, , Table 4, Rows 6 
& 7), the NGS 72 determination is in significantly better agreement with FAG 73 
than AMS 67, the area covered in these two cases being substantially equivalent 
to the United States. However, the latter solution was based on comparisons at 
a greater number of points (989 in contrast to 314). In addition, the number 
of astro-geodetic stations used in preparing both solutions is essentially the 
same. It is judged that the increased value of for AMS 67 comparisons over 
that for NGS 72 solutions is due at least in part to the use of the former in com- 
parisons at locations away from astro-geodetic stations and thereby introducing 
interpolation error. A plot of for the solution given in Row 1 of Table 4 is 
illustrated in Figure 11. 

Further to the discussion in Section 3.2, the quality of agreement in non- 
coastal areas where adequate astro- geodetic data is available, is on par with the 
results obtained for Australia (Mather, 1972b, p. 25). The reasons for the dis- 
crepancies are self-evident if Figures 1 and 14 are overlaid on Figure 11. They 
can be summarized as follows: 

e Astro-geodetic data north of parallel 60°N is too widely spaced for 
reliable determinations. 

e Gravity data off the west coast of North America is inadequate for re- 
liable gravimetric computations. 

As mentioned in Section 3.2, the large discrepancy in the Colorado region is 
due in part to the non-Stokesian effect (Fig. 12), the incompatibility of the 
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quantities being compared and the inadequate density of astro-geodetic sta- 
tions in the area. 

In view of the doubts that existed about the possibility of computing a re- 
liable astro-geodetic solution in northern Canada and Alaska because of the 
paucity of data, it was decided to investigate the orientation of NAD 27 using, 
as input: 

• astro-geodetic deflections only, north of parallel 48“N; together with 

• both undulations and deflections of the vertical south of parallel 48°N. 
Tests on this type of data distribution were carried out on the Australian Geodetic 
Datum where uniform data coverage was available (Mather, 1972c). These tests 
confirmed that stable geocentric orientation parameters could be determined 
using this type of data distribution. Solutions of this t 3 q)e are classified as 
Class C in Table 4. 

4,2.4 Solutions From Astro- Geodetic Deflections Only 

A third type of solution classified as B in Table 4 was also made in re- 
sponse to the irregular distribution of astro-geodetic data as a function of lati- 
tude (see Fig. 14) and the likelihood of systematic en'ors in any astro-geodetic 
tmdulation representation deduced from such information. These were effected 
by comparing deflections of the vertical alone against equivalent gravimetric 
values deduced from Equation 14. The resulting solutions are shown as Class B 


in Table 4 (Rows 5 and 8). 
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The following conclusions can be drawn frotn these results: 

* Comparisons over an area approximating to the United States (Row 8) 
have smaller rms residuals than those made over the entire NAD 
area. 

• The magnitude of both (Column 6) and (Column 9) are about 
2-1/2 arc sec. This is due to three factors: 

a. An inner zone of radius 10 km around the point of computation has 
been excluded when evaluating FAG 73. 

b. Non-Stokesian contributions have been excluded. 

c. Errors in excess of one arc sec can be expected in the geodetic 
coordinates of points in the data set used north of parallel 50°N 
(e.g. , Lilly, 1960). 

In fact, the values of and for the United States (Table 4, Row 8) are not 
significantly different from values obtained during similar computations in 
Australia and where the subsequent inclusion of inner zone surveys halved the 
values of and (Mather, 1970a, p. 72). 

Another observations of significance in the case of these Class B solutions 
based on the comparisons of deflections alone, concerned the changes of deter- 
minii^ AN from solutions of this lype. It has been observed in the case of 
limited extents like Australia (1-1/2% of the Earth's total surface area) that AN 
was effectively indeterminate (Ibid, p, 72, Type 2 Solutions). This was also 
observed in the case where determinations were based on comparisons confined 
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to a similar area in Nortli America (Table 4, Row 8). When the area was ex- 
tended to the entire NAD (Table 4, Row 5), a relatively realistic estimate of 
AN was obtained. This particular solution is therefore totally independent of 
any astro-geodetic und ulation solutions. It is interesting to look at the Cartesian 
components of the geocentric orientation vector (Table 4, Columns 14 through 16) 
which are in good agreement with the average satellite solution (Table 5, Row 9) 
for values of AX 2 and AX^ but not AX^ , thus apparently pointing to some degree 
of instability in the determination at right angles to both the rotation axis and 
the meridian at Meades Ranch. But also see the discussion associated with 
Table 6 in Section 4.3.2. An important corollary is that the larger the surface 
area being compax’ed, the more reliable the determination of AN from the com- 
parisons of deflections alone. 

4. 2. 5 Non-Stokesian Effects 

The effect of non-Stokesian terms is evaluated through Equation 11. These 
effects were found to have a negligible influence on determinations of the geocen- 
tric orientation parameters (see Table 4, Rows 1 & 2). The reason for this 
becomes apparent on studying Figure 12. While significant magnitudes of the 
non-Stokesian contribution to the height anomaly occur in the western part of the 
NAD, ranging from 1-3 m in amplitude, tliese are largely of short wavelength, 
the sign of the contribution being controlled by the following factors; 

« the sign of the free air geoid height Nf (Fig. 10) winch is largely 
negative in the region; and 
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• the contribution of [ sin (dh/dr) + hp - h] varying in sign as the point 
of computation moves across the Roclcy Mountains from west to east. 
Consequently there is no systematic build-up in the magnitude and wavelength 
of the non-Stokesian contribution in the western region of NAD tc bring about 
changes of significance in the geocentric orientation parameters for NAD. The 
corrections obtained are of the order of 30 milliarc sec in and 3 cm in AN 
on the basis of the present calculations. It should be emphasized that Equation 
11 is an approximation. A more complete version is given in (Mather, 1974, 
p. 100), but for reasons given in Section 1.3, it was not warranted to persevere 
with a more complex version of the expression for N,. in viev; of the excessive 
computer requirements. However, tliere is a possibility that the values of the 
non-Stokesian effect illustrated in Figure 12 may well underestimate N^, espe- 
cially in areas of rugged topography where the elevation data available in com- 
puter compatible form at the present time is not sensitive enough to model the 
fine structure of the topography. More research is needed in tliis area on a 
regional basis as point computations by themselves give no real information on 
the significance of the non-Stokesian effects in practical high precision geodesy. 

4.2.6 Conclusions 

The following conclusions can therefore be drawn from results summarized 


in Table 4. 
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i. The fi'ee air geoicl alone is adequate for determining geocentric orien- 
tation parameters for the NAD region to ±0. 05 arc sec or its equivalent 
in each coordinate and non-Stokesian terms need not be considered. 

ii. The data available at present gives comparisons between FAG 73 and 
astro-geodetic imdulations at astro-geodetic stations below parallel 
48° N wliicli an rms residual (Table 4, Colmnii 12) of ±l.Sm (Table 
4, Row 6) compared with ±1.6m for the 1971 determination for Australia 
(Mather, 1972b, p. 23), This value drops to ±1.0 m (Table 4, Row 9) 

if comparisons were restoicted to astro-geodetic stations below parallel 
40° N. The slightly larger figure than obtained in the Austx'alian study 
is attributed to the increased ruggedness of the topography on NAD and 
the resulting (correlated) incompatibility between the data types 
compared. 

iii. The difference between the geocentric orientation parameters obtained 
for the entire datum (Table 4, Row 1) and the area south of parallel 
48° N (Table 4, Row 7) illustrates the necessity for basing determina- 
tions of the geocentric orientation parameters on comparisons extending 
over the entire datum. The sigixificance of the results obtained in this 
study are, however, not entirely clear-cut in view of tJie imcertain pre- 
cision of the geodetic data in nortliern latitudes wiien taken as represen- 


tative of NAD as a whole. 
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iv. Geocentric orientation parameters detexnnined from deflection com- 
parisons alone {Table 4, Rows 5 & 8) were foiuid to give reasonable 
values of AN when based on comparisons extending over the entire 
datum, though the AN value was for all practical purposes, indetermin- 
ate when the comparisons were restricted to the United States alone. 
Neveitheless, there is a relative indeterminancy in such solutions at 
right angles to the rotation axis and the plane of the meridian at Meades 
Ranch. 

4.3 Comparison of Astro-gravimetric and Satellite Determinations of the 
Geocentric Orientation Vector 

4. 3, 1 Inter-Relation of Geocentric Oinentation Parameters 

Several satellite solutions are available for defining the displacement of the 
origin of NAD from the geocentre. Some of these are summarized in Table 5, 
along with the average astro-gravimetric determination from Table 4 for ease 
of comparison. The satellite determinations are obtained by comparing the co- 
ordinates of tracking stations as determined in global geodetic networks with 
equivalent values on the regional geodetic network using a seven parameter fit 
to allow for translation of origin (three parameters AXj), rotation of axes 
(three parameters a>j) and scale (og). Dynamic satellite solutions are generally 
linlced to a Greenwich— Conventional International Origin (CIO) system of refer- 
ence by incorporating camera (optical) data for provision of orientation. The 
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scale of the system in modern satellite solutions is usually based in part on 
laser tracldng data (e.g. , Marsh et al. , 1973; Lerch et al. , 1974) as well as 
on the value adopted for GM in the case of dynamic solutions. The scale of 
purely optical solutions is more dependent on long baselines for scale (e.g. , 
Schmid, 1974). Wliile the parameter tjg is of importance in attempting to relate 
regional geodetic datums to the geocentre, the significance of tlie rotational 
parameters coj need careful consideration when both the regional geodetic datum 
and the reference system used for the satellite solution are related to Greenwich/ 
CIO using essentially common considerations. For further discussion, see 
Section 4.3.2 and the Appendix 

The quantities Axj are simply related to the geocentric orientation param- 
eters A^j defined in Equation 26 by the formulae 

AXj = p A^ sin 0Q cos Atj sin A^ - AN cos 0^ cos A^' 

= P A^ sin 0Q sin A^ - A 17 cos A^ - AN cos 0^ sin A^ V (28) 

AX3 = -p A| cos 0Q - AN sin 0^ 

on taking into account that changes A^ and At? are equivalent to changes 50 and 
5 A in the latitude and longitude at Meades Ranch in terms of Equation 26. 

Further, a change AN in the ellipsoidal height is obtained by a radial outward 
displacement of (-AN) of the ellipsoid centre along the normal to the regional 
geodetic origin. Similarly, 

p A| = AXj sin 0 Q cos A^ + AXj sin 0 ^ sin A^ - AX 3 cos 0 ^ 

Ar? = AXj sin A^ - AXj cos A^ 

AN = -(ZiXj cos 0Q cos Aq + AXj cos 4 >q sin A^ + AX 3 sin 0 ^^] 


( 29 ) 
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4,3.2 The Effect of Scale and Axial Rotations 

All values of A|. obtained by the use of Equations 26 and 29 fi’om the input 
values of AXj for global satellite solutions are given in Table 5. The input values 
into Table 5 in the case of astro-gravimetric solutions are of course, the values 
A^j . Equivalent values on Clarke 1866 ellipsoid are obtained by the application 
of Equation 27. In assessing the equivalence of values of AXj determined from 
satellite solutions with the values of A^j determined astro-gravimetrically and 
related through Equations 28 and 29, the following points should be borne in 
mind: 

• Gravimetric determinations are insensitive to effects of zero degree 
in AN. Effects of a similar type are exlcuded from the satellite de- 
termined AX| and coiild be reflected in the scale Og. However, any 
scale difference between that prevailing on the regional geodetic datum 
and that defining the satellite solution in the case of astro-gravimetric 
determinations will be reflected in all A^j . The A Xj obtained from the 
astro-gravimetric method could thus be different from the translational 
parameters AX. obtained from satellite geodesy. 

9 This factor apart, the astro-gravimetric method is not sensitive to the 
system of reference adopted or the value of GM provided the geodetic 
origin of longitudes is the same as the astronomical origin of longitudes 
and the geodetic network is accordingly oriented in Earth space. For 
a further discussion see the Appendix. 
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• D 3 mamic satellite solutions which are related to CIO /Greenwich using 
a system similar to the concepts adopted by the Bureau International 
de I'Heure for the maintenance of the system of reference (Guinot & 
Feissel, 1969) are implicitly related to the geocentre (Earth's centre 
of mass) and therefore constitute a determination winch is geocentric 
and oriented by the CIO/Greenwich system of reference to the precision 
of the observations (±0.2 arc sec or 6 m in each component). 

It would therefore appear that it is possible to purposefully study the effect 
of introducing the rotational parameters toj in dynamic solutions when comparing 
the coordinates of tracking stations as defined on the global satellite system with 
those on the regional geodetic system. Two non-geometric satellite determina- 
tions are catalogued in Table 6, giving parameters deduced for NAD 27. The 
values of AXj for GEM 6 used in Table 6 (Lerch et al. , 1974, p. 84) are slightly 
different from those in Table 5 which were also reported elsewhere in the same 
publication (Ibid, p. 88). The other solution considered is GSFC 73 (Marsh etal. , 
1973, p. 52). In both cases, cj, and cj, are zero to within the precision of de- 
termination, in keeping with the expectations based on the discussion in the 
Appendix. In both cases, however, ujj is significantly larger than zero, being 
about one arc sec to the order of precision of the determination. This is equiva- 
lent to the occurrence of a one second discrepancy between the origin of longi- 
tudes (XjXj) plane of the satellite solution and that for NAD 27. As discussed 
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in the Appendix, this is not due to any characteristic at the Meades Ranch origin, 
but rather a function of the network of tracldng stations as a whole and represents 
a relative orientation discrepancy between the satellite and geodetic solutions 
about the X 3 axis. Such a discrepancy would only occur in astro-gravimetric 
determinations if there were an error in the orientation of the geodetic network 
on NAD 27 in relation to the CIO /Greenwich system of reference implied in the 
astronomical determinations used for computing astro-geodetic deflections of 
the vertical. If this were the case, an additional effect will be introduced, af- 
fecting primarily the quantity A 77 in Tables 4 and 5, 


The inter- comparisons of rotated satellite solutions and non-rotated astro- 
gravimetric solutions could be related as follows. The effect of counter- 
clockwise rotations about the Xj axes would produce changes 5AXj in the 
geocentric orientation parameters AXj related to the geocentric Cartesian co- 
ordinate system Xj, given by (e.g. , Mather, 1973, p. 197) 

6AXj = 

where 


f 0 if i = j or j = k or i = k 


^iik ~ r 1 subscripts occur in the order 123123. 


ijk 

I 

.-1 if subscripts occur in the order 1321321, 
for small rotations coj, Xj^ being the geocentric Cartesian coordinates of 
Meades Ranch (values on the regional geodetic datum for the coordinates of 


Meades Ranch would be adequate for practical evaluations). 
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Both GSFC 73 and GEM 6 have significant scale differences in relation to 
the geodetic network as defined on NAD 27. If the satellite solution has a scale 
wliich is greater than that of the regional geodetic network by Og parts per mil- 
lion (ppm), the primary effect is on the geocentric orientation parameter AN. 

The difference between AN^ ^ y as obtained from an astro- gravimetric determi- 
nation and ANg^^ from a satellite solution on this account is given by the relation 

^^GRAV “ 

where is the mean radius of curvature of the reference ellipsoid at Meades 
Ranch. On combining Equations 30 and 31 with Equation 28, for all practical 
purposes. 


AXi 


‘gray 


- AXi Wi Xl,,. - On X; 


‘SAT 


ijk '"j "^io 


(32) 


where AXj would refer to a set of geocentric orientation parameters on a 
•gray 

geocentric Cartesian coordinate system based on a solution involving no rota- 
tions and allowance for scale, while the value AX: would have been obtained 

‘sat 

from a seven parameter fit between a global solution and the regional geodetic 
network. 


The application of Equation 32 to GSFC 73 and GEM 6 is illustrated in 
Table 6 where comparisons are shown in relation to the average astro- 
gravimetric determination based on FAG 73. It can be seen that the largish (in 
excess of 2 ppm) discrepancies wMch exist between the geocentric orientation 
parameters AXj as determined 
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a. between the astro -gravimetric solution on the one hand and one of either 
GSFC 73 or GEM 6 on the other; and 

b, between GSFC 73 and GEM 6 

decrease by at least a factor of 2 in all cases when allowance is made for the 
rotations cjj and the effect of scale og , 

It would therefore appear that the goal of determining geocentric orientation 
parameters for the NAD 27 with the same resolution as satellite methods using 
astro- gravimetric techniques has largely been achieved. The precision sought 
of ±0.2 arc sec or its equivalent in each parameter cannot be said to have been 
achieved with confidence. At face value, the most strongly determined param- 
eter is AN. It should be aclmowledged however, that the values quoted in Table 
4 may still have a bias of up to one part per million (±6m) because the gravi- 
metric method as used in the present determination is insensitive to scale on a 
global basis. 

On the other hand, the satellite solutions are obtained by substantially in- 
dependent methods. The main dependence of the gravimetric determination on 
satellite techniques would be values for the harmonic coefficients in the GEM 4 
model with nodal points in excess of 60° apart. Any eri’ors in these values 
would influence the geocentric orientation parameters obtained by the astro- 
gravimetric method as used in tliis study. It should be noted that values ob- 


tained from the comparisons of FAG 73 with deflections of the vertical alone 



50 


(Table 4, Row 5) over the entire NAD are only in marginal disagreement (less 
than Ippm) with the eqtiivalent GEM 6 in Table 6. 

The main conclusion that can be drawn on the basis of the results in Table 6 
is that any residual distortion in the network of satellite tracldng stations on NAD 
in the solutions GSFC 73 and GEM 6 must be less than the discrepancies between 
the equivalent AX. 's as obtained for each satellite solution and the astro- 
gravimetric determination— i.e. , at the level of around Ippm. This conclusion 
is subject to the qualification given in the second sentence of the previous para- 
graph, It would certainly be of interest to repeat these tests using an astro- 
geodetic network based on the new North American Horizontal Datum (NAHD). 

One matter which remains to be dealt with is the possibility of strengthen- 
ing the values of the low degree harmonic coefficients defining the GEM 4 Earth 
gravity model. For example, information of interest is contained in Table 4, 
Rows 5 & 8. The solutions listed in these rows are independent of any errors 
in computing astro-geodetic undulations. It can be seen that the value of for 
the entire datum is smaller than that for the area below parallel 48' N by about 
0.5 arc sec. The value of A 77 , however, is almost unaffected. Two other pairs 
of solutions covering the same two areas are also listed in this table— (Rows 1 
& 7) and (Rows 3 & 6). The results are summarized on a differential basis in 


Table 7. 
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The differences between values for A| , At? and AN obtained from solutions 
extending over the entire NAD (an area subtending an angle of approximately 60“ 
at the geocentre) and those confined to the continental area below 48“ N (subtend- 
ing an angle of approximately 30“ at tlie geo centre) could be attributed to two 
possible causes provided there were no substantial systematic errors in the 
global gravity standardization networks— errors of about 1-2 mGal would be nec- 
essary to produce an effect of around 0. 5 arc sec in the orientation parameters. 
They are the following: 

• systematic errors in the geodetic network on NAD 27 wMch were vari- 
able functions of distance from Meades Ranch; and/or 

• errors in the harmonic coefficients in the GEM 4 model in the range 
3 < n < 6. 

The results in Table 7 appear to indicate that the net effect of Urjse errors 
is, for all practical purposes, zero in the east-west direction while that in the 
north-south direction is about twice the noise level of the astro-gravimetric 
determinations in the case of solutions independent of astro-geodetic undulation 
determinations (Row 1). The equations relating the changes listed in Table 7 to 
corrections in the harmonic coefficients are given in (Mather, 1970b, pp. 98-99). 
However, in view of the uncertainty associated with geodetic coordinates in the 
northern part of the NAD, a conclusive gravity model improvement is unlikely 
to be achieved using this method on the data available at present. The matter 
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should be reviewed after geodetic coordinates for all astro-geodetie stations 
used in this investigation are re-establish ed on NAHD, 

Meaningful results would only be obtained if such a determination were 
studied in conjunction with similar investigations for an equivalent area in other 
parts of the world. An investigation of this type is currently being undertaken 
for the Australian and Indonesian regions as a whole in the time frame of the 
next half decade. The combined analysis of such results with those obtained 
for the North American Region should provide confirmation of the values adopted 
for the low degree harmonic coefficients in the representation of the Earth's 
gravity field, independent of satellite values, up to degree 6. 

5. CONCLUSION 

The primary purpose of this investigation was to study the feasibility of 
establishing a world geodetic system from the comparisons of gravimetric and 
astro-geodetic determinations of the separation vector. This could be achieved 
on a step-by-step basis, considering a single regional datum at a time. The 
a priori modt 1 used for the distant zone effects of the free air anomaly field 
in the present calculations was GEM 4 any any results obtained in this pilot 
study will be biased by those errors in harmonics of degree less than four. 
Errors in the values adopted for harmonic coefficients between degrees four and 
six could, in theory, be separated and identified if a sufficient number of geode- 
tic daturas of the same extent were simultaneously investigated. The results 
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obtained in this connection from the present study have an added uncertainty due 
to widely reported distortions in the geodetic network on NAD 27 in the northern 
half of the region. 

Two sets of independent criteria are available for assessing the precision 
of the two gravimetric determinations FAG 73 and HTAN 73 (see Table 4) used 
in this investigation. 

• Firstly, the gravimetric midulations can be compared with astro-geodetic 
undulations after translational correction for the geocentric orientation 
vector for NAD 27. The root mean square (rms) residual of compari- 
sons in the normal displacement (Table 4, Column 12) are found to vary 
from ±1.0 m for the region south of parallel 40“Nto ±2. 9m for the en- 
tire datum, when defined at astro-geodetic stations representing cor- 
ners of a one degree equi-angular grid wherever such data were avail- 
able. The rms residual is slightly inferior in the case of comparisons 
where the astro-geodetic data were obtained by interpolation. 

9) The second independent criterion is the comparison of geocentric orien- 
tation as obtained from satellite geodesy with those computed by the 
astro-gravimetric method. In the former teclmique, the geocentric 
orientation parameters are obtained as three components (AXj) on a 
geocentric Cartesian coordinate system Xj as part of a seven param- 
eter transformation between coordinates on the regional geodetic sys- 
tem and the global satellite system. On allowing for the fact that any 
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scale errors in the geodetic network and any tendency for the regional geodetic 
system to be incorrectly oriented in relation to CIO/Greenwich through azimuth 
control, as discussed in the Appendix, are absorbed in astro-gravimetrically 
determined geocentric orientation parameters (A^, At?, AN) in the local Laplacian 
triad at the regional Meades Ranch origin, it is shown in Section 4. 3 that the 
average value obtained astro-gravimetrically agrees with eacn of the global satel- 
lite solutions GSFC 73 and GEM 6 better than the two satellite solutions agree 
with each other, as detailed in Table 6. 

These tests provide bounds for 

a. long wavelength (nodal points > 6000 km apart) errors in GEM 4; and/or 

b. systematic error propagation in the present geodetic network on NAD 27. 
The discrepancy between "equivalent” satellite solutions and the average astro- 
gravimetric determination listed in Table 6 also provides an upper limit for any 
zero degree term in the height anomaly. On the basis of the present study, this 
does not exceed the precision of the satellite solutions. On the assumption that 
there is no scale error in the geodetic network on NAD 27 and assuming that no 
systematic orientation error exists in relation to CIO /Greenwich as implied in 
both GSFC 73 and GEM 6, the recommended set of geocentric orientation param- 
eters for NAD 1927 on the basis of the astro-gravimetric determinations is 

A^ = 0.1 arc sec; At? = 1.0 arc sec; AN = -34 m in Laplacian Triad at Meades 
Ranch or AX, = -7m; AXj = 159 m; AX^ = 169 m on a geocentric Cartesian co- 
ordinate system. 
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Corrections which make allowance for the assumptions mentioned above on 
the basis of satellite solutions are given in Table 6. The significant parameters 
are a possible allowance for scale inconsistencies between that provided by the 
horizontal geodetic network on NAD 27 and the current values of GM (expected 
discrepancy of 1-2 ppm) and the allowance for a rotation between coordinates on 
the satellite system and the horizontal geodetic network of about 1 arc sec around 
the Xj axis. The only effect in excess of Ippm as a consequence of such con- 
siderations is a change of about -24 m in AXj and +0. 8 arc sec in At? due to the 
rotation of 1 arc sec about the X3 axis. 

The use of FAG 73 in lieu of HTAN 73 (solution for the height anomaly using 
Equation 9) is not found to make any non-negligible contribution to the geocentric 
orientation parameters for the North American Datum. This does not rule out 
the possibility that these effects could still be of significance when orienting a 
datum which included the Himalayas, The extent of the datum used has a signif- 
icant effect on the results (see Table 7) . This is tentatively attributed in the 
present case to the effects of a, and b. above. A redetermination after the es- 
tablishment of NAHD might possibly reduce the influence of b, While the use of 
the astro-gravimetric technique may only be of limited interest in future plans 
for establishing a world geodetic system, it still provides useful information for 
the study of slowly varying systematic effects and the scale of a regional geodetic 
system, while at the same time, providing a means for determiningthe geocentric 
orientation vector for regional geodetic datums independent of scale as introduced 


through satellite systems. 
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8. APPENDIX - THE ROLE OF TRANSLATIONS AND ROTATIONS IN 
ASTRO-GRAVIMETRIC DETERMINATIONS OF THE GEOCENTRIC 
ORIENTATION VECTOR 

The discussion in Section 4 shows that some ambiguity is experienced in 


practice when comparing determinations of geocentric orientation parameters 
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{datum shifts) as obtained from solutions which specifically exclude rotations 
(astro-gravimetric solutions) with those which call for seven parameter trans- 
formations, The background behind the exclusion of rotational concepts in astro- 
gravimetric methods is the following. 

The orientation of the reference ellipsoid used in gravimetric determinations 
and its location in Earth space are implicit on the basis of conditions postulated 
in the derivation of solutions to the geodetic bomidary value problem (e.g. , 
Mather, 1970b, p. 84). It can be shown that in a careful derivation of a solution 
of the geodetic boundary value problem the reference ellipsoid is located in 
Earth space so that its centre is at the centre of mass of the Earth (or, more 
specifically, the centre of mass of the solid Earth and oceans which is within 
+5 cm of the geocentre). This is due to the fact that no first degree harmonic 
can be permitted in the representation of the distrubing potential. The minor 
axis of the reference ellipsoid is placed in coincidence with the axis of greatest 
moment of inertia of the Earth if C 21 =821 = O, in the spherical harmonic ex- 
pansion of the geopotential. Actual numerical solutions for these values (e.g. , 
Lerch et al. , 19'^2) appear to indicate that their magnitude is about three orders 
smaller than other coefficients (excluding C,^) thus indicating that the effect of 
uncertainty in this regard is at the 5 cm level in solutions of the geodetic bomidary 
value problem. It can be concluded that the gravimetric solution refers to the 
CIO /Greenwich/Geocentre system of reference without ambiguity in practical 
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determinations. This system of reference can be fully represented in Earth 
space by 9 parameters; 

• The coordinates Xgj of the geocentre on some three dimensional 
Cartesian coordinate system. 

• The direction cosines 2jj of the line joining the geocentre to the CIO 
pole. 

• The direction cosines of a line in the meridian of reference for 
longitude (Greenwich) which is perpendicular to this "axis of rotation." 

If the Xj ax’S system was geocentric, the following conditions hold: 

Xgj = 0 (A-i) 

3 

Y, fifj = j = >’2 (A-2) 

i = I 

3 

E 'il 'i2 = “ 

i = 1 

It therefore follows that three parameters (e.g, , j , G, j , G, t ) have to be 
arbitrarily defined to specify a consistent system in Earth space. 

Any point P in Earth space is completely defined in relation to such a sys- 
tem of reference if its three dimensional Cartesian coordinates (Xjo) were 
Imown. Conversely, given the three coordinates (Xjy) which are arbitrarily 
assigned to the point P in Earth space, three degrees of freedom exist when 
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fixing the (unlaiown) location of the system of reference in Earth space. In the 
case where the differences 

= (A-4) 

are small (of the order of a few hundred meters), the options available can be 
visualized as follows. Consider the point Q in Earth space whose coordinates 
with x'espectto the Xj system defined by Equations A-1 to A-3 are Xj^. It can 
be seen that the assignation of the coordinates Xj^ to P could be completely de- 
scribed in Earth space by placing the origin of the rectangular Cartesian coordi- 
nate system X| at a point whose coordinates on the Xj system are AXj without 
rotation of axes . In tills concept, the quantities defined by Equations A-2 and 
A-3 have been held unaltered by choice when exercising discretionary assigna- 
tion of the three arbitrarily defined parameters. 

Needless to say, the assignation of coordinates Xj^^ to P could also have 
been achieved by rotating the axes while maintaining Equation A-1. Such a con- 
sideration is unnecessarily airJjiguous in complex geodetic inter-relations. It 
is therefore valid to consider the assignation of coordinates at a single point on 
a geodetic datum to be equivalent to a translation of origin without rotation. 

The concept of rotation is introduced to take into account systematic trends 
in a geodetic network. A classic example where a set of rotations may be of 
relevance is in instances where the geodetic network is incorrectly oriented in 
relation to CIO /Greenwich. An example of this tsqje would be a geodetic network 
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where no CIO /Greenwich related Laplace azimuth were available to orientate 
the net at the origin. Consequently, the question of rotations should not arise 
in the case of error free observations provided the Laplace azimuths introduced 
into the network were all properly linked to CIO /Greenwich using appropriate 
inertial reference and timing systems. In practice, rotations of the order of 
magnitude of errors in the geodetic network should be expected when orienting 
geodetic datums to CIO/Geocentre/Greenwich. Comparison of two such net- 
works will define rotations wliich are also a measure of the a verage distortion 
of each network. 

It must be concluded that in comparing solutions aspiring to a precision of 
Ippm, only those rotations whose magnitude exceeds 0.2 arc sec can be inter- 
preted as being due to factors other than systems uncertainty. In the case of 
solutions for the NAD 27, it is only cjj which is geodetically significant. Its 
nature appears to indicate an orientation of the North American geodetic network 
which is not consistent with CIO/Greenwich by about one arc sec as discussed 


in Section 5. 
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Table 1 

Area Sub-Divisions Used In Computations 

Ranc,e of ^ ji'ree Air Anomaly Data Type Data Sources 

(Degrees) ^ 

4/ <0,1 Circular Ring Representation Estimated from 0. 05° x 0. 05° 

area means 

0.1^’/'< 0.5 0.05" X 0. 05° area means Interpolated fi'om 0. 1° x 0. 1° 

area means 

0.5<i^<1.5 0.1° x0.1° area means Prom data banlts and 

interpolation 

1.5<i/i'<7.5 0.5° x0.5° area means From 0. 1° x 0. 1° area means 

7.0<\p<22.5 1° X 1° area means i. From 0. 5° x 0. 5° area 

means (Regions in Fig. 1) 

ii. Talwani et al. , 1972 or 
ACIC 1971 [if no i.] 

iii. GEM 4 [if no i. or ii.l 

i//>22.5 5° X 5° area means i. From 1° x 1° area means 

(Regions 1-81) 
ii. From GEM 4 [if no i.j 
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Table 2 

Free Air Anomalies - North America 

Error of Representation E{Ag}j^^° for Various Equi -Angular Squares of Size m° 


Square 

Size 

(m°) 

±mGal 

Number of / ^ 
Readings 1 

Njj Remarks 

0,5 

15.9 

113,881 

For North America (Equation 22) 

1 

17.8 

113,881 

-do- 

1 

16.5 

113,881 

For zero elevation by linear re- 
gression analysis 

1 

34.6 

4,812 

For elevations greater than 
2000m 

5 

24.8 

113,881 

For North America (Equation 22) 
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Table 3 

Free Air Anomalies - North America 
Error of Representation E{Ag}j ^ 2 ° for Half Degree Equi-angular 
Squares ds a Function of Elevation 


Elevation (h) 
(m) 

(imGal) 

Number of 
Squares 

-5000 < h < -4000 

9,5 

2 

-4000 < h < -3000 

14.9 

7 

-3000 < h < -2000 

19.1 

20 

-2000 < h < -1000 

24,8 

25 

-1000 < h < 0 

11.5 

66 

0 < h < 1000 

9.6 

7137 

1000 < h < 2000 

21.1 

815 

2000 < h < 3000 

24.1 

230 

3000 < h < 4000 

24.5 

114 
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Table 4 

North Amex’ican Datum 1927 

Geocentric Orientation Parameters From the Comparison of 
Gravimetric & Astro-Geodetic Data 


I 


L^ipliiclnti Triad ;it Moudes Uanch 


Gr:iv. 

Solution 


Astro- ^ 

Geodetic a 

Solution ^ ^0 


it 


Gcoccnlrtc 

Cartofllnn 


Liititudc 
llanRC {“> 


Longitude 

Rango (’) 


1 FAG ?a 

2 HTAN 7:i 

3 FAG 73 
■i FAC 73 
fi FAG 73 
f> FAC 73 

7 FAG 73 

8 FAG 73 
D FAG 73 

to FAG 73 
it FAG 73 
12 FAG 73 


AMS 07 
AMS 07 


NCS72&NI-:Y a -7.1 
AMS f*7 A -7. 1 

nEFI.FCTIONS n -7.2 
NGS 72 C '7.2 

COMPOSITE -7.7 

COMPOSITE -7.1 

AVKItAOE -7.1 


"At 

Nj 


"A, 

N„ 

AN 

"an 

N., 

AX, 

AX, 


Min. 

Miix. 

Min. 

Muk. 

(tSCC) 


(sect 1 

(fseci 


(in) 

rm) 


(m) 

tm) 

(ml 





fi 

7 

« 

9 

10 

11 

12 

13 

14 

15 

10 

17 

IH 

19 

20 

2.f> 

799 

1.2 

2. 8 

73-1 

1«(.0) 

4.0 

2529 

-13 

199 

177 

15 

69 

-150 

-90 

2.(5 

799 

1.2 

2.8 

731 

Ki(.ll) 

4,r. 

252(5 

-12 

U59 

177 

15 

69 

-150 

-60 

2.9 

791) 

1.0 

2.7 

734 

KM-n 

2,9 

.394 

-9 

139 

192 

29 

63 

-140 

-60 

2.(i 

799 

1.0 

2.7 

734 


1.8 

319 

-7 

154 

ir»9 

20 

4K 

-124 

-70 

2.(5 

799 

0.8 

2.7 

734 

7(-li) 

- 

- 

_2 

157 

178 

15 

99 

-iriO 

“60 

2. ] 

370 

1,0 

2.-1 

300 

Hi (-7) 

1.8 

314 

-7 

153 

1.58 

29 

48 

-124 

-70 

2.1 

370 

1.1 

O f", 

f>0G 

1«(2) 

2.C. 

989 

-12 

154 

ir»9 

29 

48 

-124 

-70 

2,2 

r.79 

0,9 

O Cj 

909 

IG(.7)< 

- 

- 

-4 

MH 

inc 

29 

4H 

-124 

-70 

2.(5 

799 

1.9 

2.8 

734 

Hi (.HI 

1.0 

203 

-0 

ir>r. 

190 

29 

40 

-124 

-70 



1.0 



ll(.Ii) 



-7 

192 

178 

15 

69 

-150 

-60 



l.O 



ll(.l) 



-7 

154 

tni 

29 

48 

-124 

-70 



1.0 



13 1 



-7 

159 

199 






Key to Tatilo •! 

Slolution Types 

FAG 73 1973 Free Air Gcolil 

HTAN 73 - 1973 HelKht Anomalies 

AMS 07 Army Map Service Astro-GerMletlc Solution 
(Fischer et nl. , 19R7) 

NGS 72 t*.S, National Gomiutle Survev Astfn-lieodetle 

LevoIUng iKiev. 1972) 

NEV Canadian Astro- Geodetic LcvcUitik (NEY t9;>2i 


Cliiss 

A 
B 
C 


Undulation Compartsons Only 
Ucnection Comparlsone Only 

An at B plus undulation comparisons south of Latitude in 
Column 1» 


o rni» Hesiduiil 

N Numl)or of Comparisons Made 

• Value Obtainwi 019.91 replaced by best Value from Rov.' C 
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Table 5 

Geocentric Orientation Parameters for the North American Datum 1927 


R 



In Laplacian 

Geocentric 

e 



Triad At 


Cartesian 

f 

Solution Name 

Source 

Meades Ranch 

Components 

N 




Atj 

an 

AXj 

AX 

AX 3 

0 



(sec) 

(sec) 

(m) 

(m) 

(m) 

(m) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

GEM 4 

Lerch et al. , 1972 

-7.6 

1.5 

0 

-24 

153 

181 

2 

GEM 6 

Lerch et al. , 1974 

-7.8 

1.5 

-2 

-22 

155 

187 

3 

GSFC 73 

Marsh etal., 1973 

-7.7 

2.2 

6 

-43 

162 

179 

4 

DOD WGS 72 

Seppelin, 1974 

-7.5 

1.5 

6 

-22 

157 

176 

5 

£E III 

Gaposchldn, 1973 

-7.4 

1.7 

3 

-31 

154 

176 

6 

WN 14 

Mueller, 1974 

-7.4 

1.1 

2 

-11 

148 

175 

7 

NGS 

SchmJd, 1974 

-7.7 

1.8 

18 

-32 

171 

173 

8 

Average 

Non-Geometric 

-7.6 

1.7 

3 

-28 

156 

179 

9 

Average 

Geometric 

-7.5 

1.5 

9 

-22 

159 

174 

10 

Average 

Astro-Gravimetric 

-7.4 

1.0 

14 

-7 

159 

169 








Table 6 


Comparison Between Satellite and Astro-Gravimetric Determinations of 
Geocentric Orientation Parameters for North American Datum 1927 After 
Allowance for Rotation and Seale 


Description 


GSFC 73 


GKM 6 


Satellite Value 
Scale (o^ ) ppm 

r 

(uj| >arc sec 

a Rolatioti'^itoj )arc sec 
r 

|(cj^ )arc sec 
K(|uiva!ent Satellite 
Averaijo Astro-Gravimetric 
Astro-srav. Minus Satellite 


0.»i 


SB. 



Geocentric 

In Laplacian 


Geocentric 

In Laplacian 

& 

V 

Cartesian 

Triad at 


V 

Cartesian 

Triad at 



a 

1 

Components 

Meades Ranch 

a 

1 

Components 

Meades Ranch 


u 

c 

AX, 

AX, 

AX, 

A? 

At? 

m 

u 

e 

AX, 

AXj 

AX, 


At? 

AN 

S 


im) 

(m) 

(m) 

(sec) 

(sec) 

(m) 


(m) 

(m) 

(m) 

(sec) 

(sec) 

(m) 

t?. 


-43 

162 

179 

-7.7 

2.2 

6 


-24 

151 

187 

-7.7 

1.5 

-5 


.9 

1 

4 

-4 




1.7 

1 

8 

-7 





.05 

- 

-1 

-1 




0.2 

- 

-4 

-5 





^ 2 

4 

- 

1 




-0.1 

_2 


-0 





.1 

26 

-4 

- 




0.8 

19 

-3 

- 






-12 

161 

175 

-7.6 

1.2 

11 


-6 

152 

175 

-7.4 

0.9 

5 



-7 

159 

169 

-7.4 

1.0 

14 


-7 

159 

169 

-7.4 

1.0 

14 



5 

-3 

-6 

0.2 

-0.2 

3 


-1 

7 

-6 

0.0 

0.1 

9 



Rms Discrepancy 


Astro-grav minus GSFC 73 = ±7.3m 


Astro-grav minus GEM G = t8, 7 m 

GSFC 73 - GEM G = ±10. 8 m (Equivalent); (±25.0 if uncorrected) 


72 


Table 7 

The Effect of Area Utiliaed in Astro-Gravimetric Solutions on the 


Geocentric Orientation Parameters for North American Datum 1927 


Astro-geodetic 
Data Type 

Class of 
Solution 

Row 

Reference 

in 

Table 4 

Change in Parameters on Reduction 
of Area from Entire Datum to Area 
South of Parallel 48° N 

Change in 
Af (sec) 

Change in 
An (sec) 

Change in 
AN (m) 

Deflections 

Only 

B 

5,8 

-0.5 

0.0 

Indeterminate 

AMS 67 

A 

1,7 

-0.7 

0.1 

-0.2 

NGS 72 & 
NEY 

A 

3,6 

-0.2 

0.0 

0.0 
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Figure 1. 
Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 
Figure 7. 


Figure 8, 


FIGURE CAPTIONS 

North America — Area Sub-Divisions for Gravity Data Processing 

North America — Distx'ibution of Available Free Air Anomaly Data — 
n = number of stations per 1® x 1° square 

North America — Representation of Free Air Anomaly Data Set (1° x 
1° means) GRS 67 — Contour Interval 20mGal 

North America — Topography From the Available 1“ x 1° Mean 
Elevations — Contour Interval 500 m 

Frequency Histogram Showing Occurrence of One Degree Equi- 
Angular Elevation Means — North America 

Free Air Anomaly and Elevation Correlation with Latitude 

North America — Frequency Histogram Showing Occurrence of E{Ag} 
for One Degree Equi-angular Squares 

E{Ag|jofor North American Datum - Free Air Anomalies Correla- 
tion Characteristics of the Available Sample with Elevation 


Figure 9. North America — Error of Representation for 1° x 1° Squares 
E{Ag} (±mGal) 
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Figure 10, 


Figure 11. 


Figure 12. 


Figure 13. 
Figure 14, 


North America — Free Air Geoid 1973 (GRS 67) — Contour Ihtervai 
2m 

North America — Discrepancy Between Astro-Geodetic (AMS 67) & 
FAG 73 Determinations — Contour Interval 5 m 

North America — Non-Stokesian Contribution to Height Anomaly 
(i//<5‘’) “ Contour Interval Im 

Relation Between RS 1967 and NAD 27 in Meridian of Meades Ranch 

North America — Distribution of Astro-Geodetic Stations 
n = Number of Stations — per 1° x 1° Square 





Fijrure 2. North America — Distribution of Available Free Air Anomaly Data 
n = number of stations per 1® x 1° square 
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Figure 3. North America — Representation of Free Air Anomaly Data Set (1° x means) 

GRS 67 Contour Interval 20 mGal 
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Figure 5 . Frequency Histogram Showing Occurrence of One Degree Equi-Angular 

Elevation Means — North Amex’ica 



Figure 6, Free Air Anomaly and Elevation Correlation with Latitude 
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I-'ifTure 9. North America — Krror of Representation for 1“ x 1° Squares E{Ag}(±mGal) 











Figure 12. North America — Non-Stokesian Contribution to Height Anomaly (i//<5“)t — 

Contour Interval Im 
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Figure 13. Relation Between RS 1967 and NAD 27 in Meridian of Meades Ranch 





